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Abstract 

We discuss the notion of integrability in quantum mechanics. Starting 
from a review of some definitions commonly used in the literature, we 
propose a different set of criteria, leading to a classification of models in 
terms of different integrability classes. We end by highlighting some of the 
expected physical properties associated to models fulfilling the proposed 
criteria. 
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1 Introduction 



Classical mechanics is a subject with a unique level of maturity: it is one of 
the most enjoyable to learn, and through the beauty and powerfulness of its 
formalism is often considered as a prototypical example of 'how things should 
ideally be done' in physics. One of the most powerful concepts in the study 
of the dynamics of classical systems is the notion of integrability (in the sense 
of Liouville, see e.g. [T]), namely that if a system with n degrees of freedom 
(i.e. with 2n-dimensional phase space) possesses n independent first integrals 
of motion in involution (i.e. Poisson-commuting) , then the system is integrable 
by quadratures. The meaning of 'integrable' here is thus transparent, namely 
that the differential equations describing the time evolution can be explicitly 
integrated using action-angle variables. The solutions of the equations of motion 
thus display periodic motion on tori in phase space, and ergodicity is absent, 
in contrast to non-integrable models which explore phase space densely in the 
course of time. Besides providing explicit solutions to the time evolution, the 
classical notion of integrability thus partitions classical models into separate 
classes of integrable and non-integrable models with manifestly different physical 
behaviour. 

It thus comes as a surprising (and insufficiently known) fact that trans- 
lating the notion of integrability to the quantum context has faced numerous 
pitfalls, and remains to this day a subject of debate. This leads to some un- 
fortunate widespread confusion, since integrability is mentioned very often in 
contemporary discussions and publications concerning among other themes in- 
and out-of-equilibrium dynamics, relaxation and thermalization of many-body 
quantum systems under current theoretical and experimental investigation. If 
quantum integrability is ill-defined, how can we thus invoke it at all? 

Questioning the precise meaning of 'quantum integrability' has been done 
on many occasions. Nearly two decades ago, in an eminently readable article, 
Weigert [2] summarized some fundamental issues and discussed the shortcom- 
ings of commonly used definitions. Delving further into the details is however 
not usually done in research articles, but rather in private discussions or pro- 
ceedings of lectures given by eminent researchers in the field, e.g. {3^. It is also 
a subject of ongoing work (see for a recent example ^). Since quantum integra- 
bility was very often mentioned and discussed by many participants during the 
StatPhys 24 conference, we found it appropriate to use the occasion offered by 
these proceedings to share a few hopefully worthwhile thoughts, observations, 
suggestions and conclusions on this important theme. 

The paper is organized as follows. We first put the problem in context, 
highlighting precisely what the problem is, and what kind of solution would 
ideally be required. In Section [3l we review many definitions commonly found 
in the literature, and collect our thoughts and comments on each of them. 
After summarizing the conclusions reached, we propose a new categorization in 
Section|4l and provide examples of where known models fit within our scheme in 
Section [5j The physical consequences of our definition are discussed in Section 
[6l which is followed by our conclusion. 
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2 Motivations 



We can begin with the simplest amd most important questions: Don't we have 
a proper definition already? Why is this question interesting and important? 
The first question will be answered in the negative in the following section. The 
second question is best answered by reflecting on the classical case: since the 
presence or absence of integrability in a classical system is associated to such 
drastic differences in physical behaviour, the lack of a proper understanding 
of the quantum equivalent inevitably means that we must be 'missing out' on 
some important properties and features; the lack of a quantum equivalent to 
the KAM theorem [3 [6j (on the stability of quasi-periodic motion in the 
presence of small perturbations) is possibly the most striking illustration of this 
point, and makes it difficult to extract hard statements on the equilibration and 
thermalization of many-body quantum systems. 

The lack of correspondence between classical and quantum integrability, 
which we will discuss further below, leads us to ask more basic questions about 
the differences between classical and quantum systems ^Sj . A first point worth re- 
membering is that quantum mechanics differs markedly from classical mechanics 
in the way it counts degrees of freedom. In quantum mechanics, discretization 
of levels means that we can comfortably work with finite-dimensional Hilbert 
spaces: spins, bound atomic levels, etc. have eigenstates which we can label with 
a discrete quantum number taking a finite set of values, and we typically say 
that the number of degrees of freedom of a quantum system is the dimension- 
ality of its Hilbert space. By contrast, in classical mechanics, we count degrees 
of freedom by specifying how many pairs of conjugate phase space variables 
are necessary to specify the configuration of a system. Each variable can take 
on a continuum of values. In any quantum-classical correspondence, we would 
thus associate the number of classical degrees of freedom to the multiplicity of 
infinities of the dimension of the Hilbert space. There thus cannot be a clas- 
sical equivalent to a quantum system with a finite-dimensional Hilbert space, 
and this already means that classical integrability is insufficient as a basis for 
defining quantum integrability in general. 

When thinking about conserved charges, the notion of Liouville integrability 
includes a specification of how many independent charges we need, namely a 
number identical to the number of degrees of freedom n (in which case it is said 
that the system possesses a complete set of charges, or is completely integrable). 
If we can provide more than n charges, the system is said to be superintegrable; 
if we can produce 2n charges, the system is maximally superintegrable (see 
e.g. [H] and references therein). One semantic pitfall is associated to the word 
'complete'. Namely, one fundamental notion in quantum mechanics is that of 
a complete set of commuting observables (CSCO), namely a set of commuting 
quantum operators whose eigenvalues are sufficient to uniquely specify a state 
in Hilbert space. In the context of integrability, the word 'complete' takes 
on a different meaning: for a fully nondegenerate system, a single operator (the 
Hamiltonian) already forms a CSCO. The cardinality of a CSCO is thus patently 
not the number of conserved charges we should be looking for in the quantum 
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case. We should of course be looking for a maximal abelian subalgebra of 
quantum operators in Hilbert space, meaning that we should be able to display 
a number of commuting operators coinciding with the dimensionality of the 
Hilbert space in order to call our set of charges 'complete'. To avoid this pitfall, 
we will thus rather talk about 'maximal' sets than 'complete' ones. 

Before going further with our discussion of quantum integrability, it is worth- 
while to follow the example of j2] and formulate a number of requirements for 
a meaningful and useful definition of this concept. Most importantly, 

1. it should be unambiguous; 

2. it should partition the set of all possible quantum models into distinct 
classes; 

3. different classes of models should display distinguishable physical behaviour. 

In addition to these, we could formulate a number of extra requirements, namely: 
a) the contact with the classical limit should be natural; b) the contact with 
integrable field theory should be natural; c) the different classes should be (algo- 
rithmically) distinguishable, i.e. it should be easily feasible to determine which 
class a model belongs to, etc. These are however less crucial criteria than the 
ones we have selected. 

3 Common definitions used in the literature 

In this section, we initiate our discussion by summarizing a number of definitions 
of quantum integrability encountered in the literature. We briefly comment each 
one. 

QI:N A system is quantum integrable QI:N if it possesses a maximal set of 
independent commuting quantum operators Qa, a = 1, dim{'H). 

Allowing for a bit of flexibility in the precise terms used, this is (at least 
in spirit) overwhelmingly the most common definition of quantum integrability 
encountered in the literature. It has the appeal of being directly related to the 
classical notion of integrability, in the sense of being essentially a word-for-word 
translation after replacing Poisson brackets with commutators. 

Definition QI:N is given the label N for a simple reason: it is too naive. 
Its fatal flaw is absolutely trivial: all quantum models associated to (limits 
of) finite-dimensional Hilbert spaces fall under the label QI:N. By the spectral 
theorem, all Hermitian Hamiltonians are readily diagonalizable; one thus obtains 
dim(?^) orthogonal state vectors \'^a) from which one can build projectors Qa = 
|^'q)(^q|, the set of which constitutes a maximal independent commuting set. 
So is every quantum system we can think of to be called integrable? Well, a 
court jester would amuse himself doing precisely this, but we have to reject this 
pathway by invoking one of the requirements we had about a proper definition, 
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namely that it should separate models into distinct classes. Definiton QI:N 
blatantly fails to do so, and is thus to be rejected as being formally useless. 

Loopholes of definition QI:N were also extensively discussed in [2]. Loop- 
hole (A) there corresponds to the flaw mentioned above. A second loophole 
mentioned is associated to the notions of 'maximal' and 'independent' quan- 
tum operators: a theorem of von Neumann [TU] states that it is possible to 
encode any number of commuting Hermitian operators into a single Hermitian 
operator Q (in other words, any operator Qq, can then be viewed as a function 
Qa = /q(Q)), so the very basic notion of the number of independent operators 
actually seems ill-defined. Going further, since we allegedly can't even prop- 
erly count the number of charges we have, we are then prevented from honestly 
declaring that a set be maximal. 

While we cannot repair the first fatal flaw of definition QI:N, let us address 
the second point in more detail. We do not dispute the validity of von Neu- 
mann's theorem, however we do not agree that it is of relevance here. First 
of all, thinking about the counting of conserved charges which can be defined 
as how one would represent them in terms of matrices in the eigenbasis, there 
is no doubt that the number required to call a set maximal coincides with the 
dimensionality of Hilbert space dim('H), since this is the number of indepen- 
dent diagonal entries. As far as counting is concerned, the notion of algebraic 
independence (that is, the set of charges does not obey any nontrivial polyno- 
mial equation) is sufficient to make it well-defined, and is already in use in the 
literature (see for example [9] and references therein). 

We thus have to look for something beyond the naive definition. Of course, 
one of the main assumptions was that in practice, we could actually diagonalize 
the Hamiltonian to obtain the charges as projectors. This suggests a more 
pragmatic definition: 

QI:ES A system is quantum integrable QI:ES if it is exactly solvable, in other 
words if we can construct its full set of eigenstates explicitly. 

While this reminds us of the action-angle variables in the classical case, the 
reader will probably agree that this washed-down definition does not take us 
very far. This definition could be further categorized according to which method 
is employed to obtain the eigenstates: Fourier transform for free theories, Bethe 
Ansatz for specific models (although the completeness of the set of Bethe eigen- 
states is not formally proven for all models), etc. We also reject this definition 
on the grounds that it does not fulfill all our criteria: the third, in particular, 
is hard to relate to. 

QI:HO A system is quantum integrable QI:HO if it can be mapped to har- 
monic oscillators. 

This is not really practical: such a mapping is hard to construct explicitly, 
even for models we know how to solve exactly. Anyway, the existence of such a 
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mapping is guaranteed by the fact that ah C* algebras are unitarily equivalent. 

A related definition is the one used in [11], where spin systems are classi- 
fied as integrable if there exists an operator mapping leading to a Hamiltonian 
expressible in terms of action operators (i.e. J^)- 

In the mathematical literature, there exists a so-called algebro-geometric 
definition of quantum complete integrability |12) . implying that a system is 
algebraically integrable {i.e. its eigenvalue problem is solvable by quadratures, 
mirroring the classical case). This allows nice things such as the application of 
the logic of Darboux transformations in the quantum context (see e.g. [13]) to 
generate nontrivial models. This definition applies naturally to quantum models 
with continuous degrees of freedom (the notion of 'complete' is then like in the 
classical case), but we find it to be too close to the classical definition to be 
sufficiently general. 

Instead of looking at exact solvability, the way out might be to look more 
precisely at the nature of the wavefunctions. This is in particular what the 
Bethe Ansatz teaches us: the scattering in one-dimensional integrable theories 
has the remarkable feature of being factorizable in two-body scatterings. This 
forms the basis for Sutherland's definition: 

QI:ND A system is quantum integrable QI:ND if the scattering it supports 
is nondijjractive 114^ . 

This definition is very appealing. Although it does not directly mention the 
existence or the number of conserved charges, one can easily construct those 
from e.g. the set of conserved momenta. It relates more or less directly to the 
classical definition, but is also well-defined, testable (at least in principle), and 
physically meaningful in the sense that nondiffractive scattering can be trans- 
lated directly into non-ergodicity. It encompasses relativistic integrable field 
theory, where this nondiffraction condition is simply the factorization condi- 
tion. 

One shortcoming is that it relies on observing the effects of scattering in the 
'asymptotic region' in the space of states (where all particles are clearly sep- 
arated), and is thus only suited for models defined directly in the continuum. 
While we find this definition the best one available, we would rather find non- 
difi^ractive scattering as a consequence of integrability, rather than as its defining 
feature. 

We can go further and try to formulate other 'physical' definitions, for ex- 
ample 

QI:ELS A system is quantum integrable if its energy level statistics is Pois- 
sonian. 

Looking at the energy level statistics of quantum models is another very 
appealing way to address the question of integrability, since it directly connects 
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with classical mechanics via Berry and Tabor's semiclassical reasoning [TS]. This 
showed that the energy level statistics of generic quantum integrable systems is 
Poissonian. Since this is based on semi-classical reasonings, only quantum sys- 
tems with continuous degrees of freedom (and thus infinite-dimensional Hilbert 
spaces) are concerned. This was however further investigated for quantum lat- 
tice models [HI [T71 [TB] such as Heisenberg, t — J and Hubbard-like models, 
explicitly verifying that for nonintegrable cases the statistics becomes Wigner's 
GOE. For Richardson-Gaudin-like models, the degree of freedom afforded by the 
large number of tunable internal parameters means that essentially any distribu- 
tion can be mimicked; however, these quickly become Poissonian upon turning 
the interaction on [TO]. We also refer the reader to the recent preprint [50] for 
recent examples and more extensive literature citations. A related definition is 

QI:LC A system is quantum integrable if it shows level crossings (i.e. does 
not show level repulsion). 

Here, the absence of level crossings in e.g. numerical solutions of finite mod- 
els is then interpreted as an observation of non-integrability (see for example 
[21] and references therein) . This definition really pertains to families of models 
with a tunable parameter, and again cannot be sufficiently general. 

Overall, we would rather view these features of energy level distributions 
used in QI:ELS and QI:LC as consequences of integrability, rather than as def- 
initions. In any case, there exist counterexamples, for example, the Richardson- 
like models mentioned above, and Haldane-Shastry type models [22l [23l [24l [25l 
[26l [27l [28] , the main conclusion there being that the energy level statistics was 
neither Poissonian nor of Wigner type. There is thus no universally acceptable 
definition of the form QI:ELS; on the other hand, QI:LC formally only makes 
sense when discussing models with tunable parameters, which cannot cover all 
cases of interest. 

3.1 Where does this leave us? 

Let us attempt a summary of the important lessons of the preceding discus- 
sion. Most importantly, the mere existence of conserved charges (which for any 
model with a finite Hilbert space is a trivial consequence of the spectral theo- 
rem) cannot be a sufficiently meaningful criterion to classify quantum models 
in different classes. It is therefore imperative to look into the actual structure 
of the conserved charges when trying to distinguish classes of quantum models. 
But the structure of quantum operators really depends on the representation 
we use to write them down. We can always choose to work in the Hamilto- 
nian eigenbasis, in which the charges are diagonal, and we can even define them 
such that they have a single matrix entry (this is the case of the projectors 
onto eigenstates). Therefore, the existence of a basis in which the charges look 
simple again does not provide us with a classification criterion, and a practical 
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definition of quantum integrabiUty cannot exist which does not identify a 'pre- 
ferred' basis in which the structure of the charges must be looked at. In view of 
our (admittedly, optional) desire to keep a semblance of quantum-classical cor- 
respondence, we could restrict e.g. to 'real space' and 'momentum space' bases, 
or more generally accept that a definition of quantum integrability is only valid 
in a specified basis. 

As we have mentioned before, the counting of the number of conserved 
charges is also a subtle point leading to various pitfalls, as our discussion about 
von Neumann's theorem illustrated. A further pitfall occurs when counting is 
mentioned together with a physically meaningful concept, that of locality: there 
is a simple contradiction in terms when talking about a 'maximal' (or the un- 
fortunately more commonly used 'complete') set of 'local' conserved charges. 
Enforcing 'locality' restricts the number of operators which can be defined to a 
number which is negligible in comparison to the required number for 'maximal- 
ness' (the dimensionality of the Hilbert space). It is thus patently impossible 
to have a maximal set of local charges, and one of the two concepts has to give 
way: either our required set is not maximal, or it contains nonlocal charges. 
This confusion probably originates from careless interpretations of the inte- 
grable quantum field theory literature, by confusing an infinite set of charges 
with a maximal one. This is also mentioned in the context of (the infinite-size 
limit of) integrable lattice models: for example, in the fundamental paper ^5] . 
the set of conserved charges for the infinite Heisenberg chain was shown to form 
a maximal abelian family, and this was equated (first paragraph of main text) to 
an ' explicit, complete set of mutually commuting operator invariants which have 
local densities...'. We take issue with this use of words, and prefer to rephrase 
this in terms of the weaker statement that there then exists an infinite number 
of conserved charges which can be written as integrals of local densities, without 
specifying that this set is 'maximal' (it is not; it only is 'maximal' with respect 
to local operators). This leaves us with the question: how many charges do we 
need? The lesson from the earlier discussion is that we can probably drop the 
'maximality' requirement and still make meaningful statements. 

With this in mind, let us now attempt to formulate a definition which takes 
care of all these issues. 

4 Alternate definition of quantum integrability 

We define a size sequence as an infinite sequence of strictly increasing integers 
{Ni,N2, N3, ...), Ni < N2 < N3 < .... To each given size Na in the size sequence, 
we associate a Hilbert space H*^^"-* obtained from tensoring Na elemental Hilbert 
spaces Hj, j — 1 , . . . , Na . We assume that each elemental Hilbert space is finite- 
dimensional, dim(Hj) = dj < 00, so dim(H^^°^) = Jlj^i "^j = d^^"^ is also 
finite. As we move up the size sequence, we assume that we simply tensor in 
additional elemental spaces, = H^^") ®j^jv1-hi '^M 

^Contrary to what might be assumed, there is no implicit restriction to one-dimensional 
systems. 



8 



Operators in T-Lj can be represented by dj x dj Hermitian matrices, which can 
be decomposed in a basis ej, i = 1, d^. Operators Q(^°) in H'-^-^ can thus be 

decomposed in the {d^^"^)'^ basis matrices e'i-- *"° = O^ie^^ (which we call the 

preferred basis), Q(^") — J^i^ i„ Qi^^l^^ e*i---*"». For a given operator Q in a 

given preferred basis e'^ "'"" , we denote the number of nonzero entries Q['^'^\^ 

as Ne(Q'^"'). For a size sequence of operators (Q^^^), Q'^"), Q^^^), ...), we 
define the density character as the nature|j of the minimal function f{Na) such 
thatNe(Q(^"))</(^a)Va. 

We can now proceed to define a corresponding size sequence of Hermitian 
operators {H'-'^'-\ H'-^^\ H'-'^^\ ...), which we interpret as Hamiltonians acting 
in their respective Hilbert space T-L'^^'^\ We assume that the Hamiltonians in 
the size sequence are built according to some simple and meaningful algorithm 
A : 'H^^"^ —7- H^-^'^\ We restrict to Hamiltonians having polynomial density 
character, in other words to Hamiltonians which can be expressed in terms of a 
finite number of sums over elemental labels. 

Being Hermitian and finite, each Hamiltonian H^^'^'> obeys the spectral the- 
orem and therefore possesses a complete set of eigenvectors. This also means 
that it automatically possesses a maximal set of conserved charges {Si^"''}, 
a = 1, d^''^''-' in involution. We can arrange the labelling of the charges in 
such a way that meaningful finite size sequences {Q[^°'\ qL^°*^\ qL^°^^\ ...), 
for all individual qL^"^ with a = 1, amax ^ can be defined. 

Definition: we will call a Hamiltonian H 0{f{N)) quantum integrable if it is 
a member of a sequence of operators having 0{f{N)) 

density character in the preferred basis, for which it is possible to define a 
sequence of sets of operators ({Q^^^)}, {Q^^^^}, {Q^^^'H, ...) such that 

1. all operators Q^^"^ in {Q'^°'} commute with each other and with their 
Hamiltonian i/'^"'; 

2. the operators in {Q*^^"-*} are algebraically independent; 

3. the cardinality C^^"^ of the set {Q'-^"^} becomes unbounded in the infinite 
size limilH; 

4. each member a — 1,...,C^^"^ of the set {Q^^"^} can be embed- 
ded within a sequence of operators {Q\^^\ Si^'', Qi^'\ ...) with 0{f{N)) 
density character in the preferred basis. 

We can thus talk about models in the linear, polynomial, quasi-polynomial, 
etc. integrability classes. For the preferred basis, while our formulation is 

■^In our definition, the terminology we use for 0{f{N)) is identical to that used for classify- 
ing algorithmic time complexity. We can thus talk about linear, sublinear, quasi-polynomial, 
suhexponential, etc. If needed, a term proportional to the unit operator should be subtracted 
when establishing the density character of an operator. 

''Notice the fact that we do not require the set of conserved charges to be maximal. 
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purposefully left generic, we are of course mainly thinking about a real-space 
based basis (in this case, the Fourier basis is also automatically compatible). 

The reader might feel that this definition falls a bit 'out of the blue'. We 
will attempt to make its meaning clearer by first specifying which classes known 
models fall into, and thereafter discussing in more physical terms how this def- 
inition should be understood. 

5 Which integrabiUty class do known models 
belong to? 

Having suggested to classify quantum models into different categories, our task 
now turns to giving explicit examples. 

Free theories 

There is little doubt that free theories should be considered integrable in some 
way, and that they should be the simplest types of integrable models available. 
Imagine thus that we have a lattice model for a set of particles, with specified 
hopping amplitudes (with no restriction on dimensionality or on these being 
local and/or long-range) and without any form of interaction. Diagonalization 
can then be trivially achieved via Fourier transformation. Conserved charges 
can be directly built from the operators representing occupation of the Fourier 
modes. If we imagine for example a size sequence in which we double the system 
size at each step and always use periodic boundary conditions, we can embed 
each such conserved charge into a sequence of constant density character. We 
thus say that free hopping Hamiltonians are constant quantum integrable in the 
Fourier basis. Looked at in the real-space basis, we can simply Fourier transform 
the conserved charges: these then become 0{N), and we can thus say that free 
hopping Hamiltonians are linear quantum integrable in the sites basis. 

EXindamental models from quantum inverse scattering 

Going beyond free theories, the most straightforward class of models to discuss 
are those which fit within the scheme of the quantum inverse scattering method 
(see [3D] and references therein), in other words which are treatable with the 
Algebraic Bethe Ansatz. This relies on the existence of i?- matrices solving the 
Yang-Baxter equation 

i?12(A,^)i?13(A, l^)i?23(A*;'^) ^ R23.{P;l')Rl3.{\y)Rl2{\y)- (1) 

The existence of monodromy matrices T{X) obeying the intertwining relation 

R^2{\,^Ji)T^{X)T2{^i)=T2{p)T^{\)R^2{\^i) (2) 

then allows to construct integrable models and their conserved charges by tak- 
ing logarithmic derivatives of the transfer matrix t(A) = Tr T(A): the set of 
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operators 

Qa = c„^lnT(A)U=4 (3) 

(where ^ is some fixed evaluation parameter and Cq. are arbitrary numerical 
constants) automatically constitutes an abelian set by virtue of ([2]). It is rea- 
sonable to assume that the set of charges so obtained is maximal, though this 
remains a conjecture. This logic allows to construct integrable quantum lattice 
models with local Hamiltonians by starting from a monodromy matrix which is 
a product of local operators, for example Heiscnberg chains [31 and more gen- 
eral models [351 [33] . Interestingly, the search for a classification of all possible 
solutions to the Yang-Baxter equation was initiated in works such as [331 (33 
but remains open to this day. 

For the explicit case of the Heisenberg spin chain, 

N 

Hxxx = JY.^,-^j+i (4) 

i=i 

(which has linear density character) the locality of the simplest conserved charges 
in the infinite system (here, locality means that the conserved charges can be 
expressed as integrals of local densities in the continuum limit) was proven in 
|36] , the structure being explicitly of the form 

Qn= ^ G^_i(ji, j„_i), 

Gn{jl,-;jn) = 0, \jl - j„ | > fl. (5) 

where the summation is over ordered subsets {ji, . . . , j„-i} of the chain and 
is a translationally invariant function. The set of local charges (which, we remind 
the reader, cannot be the maximal set of conserved charges if n/N — > where 
N is the total number of sites) was shown in [53] to be maximally abelian on 
the subspace of states where all but a finite number of spins point in the same 
direction. For finite sytems, the detailed structure of the conserved charges 
was systematically studied in a series of papers by Grabowski and Mathieu 
[371 [351 13^ , and the explicit structure was found both for periodic and open 
boundary conditions. For example, for the XXX chain, the simplest charge 
above the Hamiltonian is 

Q3=^S,-(S,+i xS,+2). (6) 

j 

A very interesting further paper by the same authors ^38] suggests to use the ex- 
istence of a single higher local conserved charge as a litmus test for integrability. 
The conjectures offered in this paper are that 1) a necessary condition for inte- 
grability is the existence of a conserved charge made of terms coupling at most 
three sites, and 2) a sufficient condition is the existence of a boost operator B 
(in Baxter's logic g^) generating higher charges according to [B, Q„] = Qn+i- 
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The authors recover many known models using this logic, but while their ideas 
are very appealing and suggestive, they unfortunately didn't produce new pre- 
viously unknown cases. 

By looking at the structure of the conserved charges above, we can immedi- 
ately see that each model obtained from a transfer matrix composed of products 
of quasi-local operators has conserved charges having linear density character 
in the sites basis. All of these are thus linear quantum integrable in the sites 
basis. Note that a similar reasoning holds for inhomogeneous models, which can 
even be used to simulate disordered systems IH] . The list of what we call linear 
quantum integrable models thus includes many famous lattice models such as 
Heisenberg chains (also of higher spin), the t — J and Hubbard models. It also 
includes variants such as restricted Bose-Hubbard with up to two bosons per 
site, or low-density limits of such lattice models yielding continuum counterparts 
like the Lieb-Liniger model 

Haldane-Shastry-type models 

An interesting generalization of the Heisenberg chain is the Haldane-Shastry 
long-range interacting model [43l |44] , 

N ^ ^ 

Hhs = X! T^^~r~^n ■ (7) 

in which Zj = e^'^'^H^ and Zj^j^ — Zj^ — Zj^ . Since it involves a double sum over 
sites of simple operator combinations, the Hamiltonian of the Haldane-Shastry 
model is now an operator with quadratic density character. The conserved 
charges of the Haldane-Shastry chain can be constructed explicitly [IS]. The 
first one above the Hamiltonian can be written 

23 - ^ E T^v'^S,, . (S,, X S,3). (8) 

Since this charge involves three space summations, it belongs to a size sequence 
with cubic density character. Each successive conserved charge includes one ad- 
ditional summation over sites. A specific conserved charge Qn thus has density 
character 0{N'"-). Some care must thus be taken to characterize the class of 
integrability this model belongs to according to our definition. The cardinality 
of the set of conserved charges becomes unbounded as — ?• oo, but a restricted 
set of charges up to some Qn^^w '^^'^ only be embedded in operator size se- 
quences of density character ©(A^""""^). To fulfill condition 3 of our definition, 
the Haldane-Shastry model must thus display a form of integrability which is 
distinct from the previous fundamental lattice models (like e.g. the Heisenberg 
chain): it cannot be a polynomial quantum integrable model. We can fit our 
definition of quantum integrability by choosing for example the density char- 
acter 0(A^'"^) ^ 0(e'" ^), which would make the Haldane-Shastry chain a 
quasi-polynomial quantum integrable model. 
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The Haldane-Shastry chain, which according to our definition thus displays 
a weaker form of integrability than e.g. the Heisenberg model, is part of a more 
general family of long-range models of Calogero-Sutherland type [351 SZl EH], 
containing many generalizations [IHl [SHI [SI]- Note that this class of models 
contains Heisenberg spin chains as a simple limit, and also more generally fits 
within the Yang-Baxter logic [S5], but the transfer matrix now gives rise to 
the generators of a Yangian, a set of non-local charges which commute with 
all local conserved charges but not among themselves. Quasi-local conserved 
charges, including the one given above, can be obtained by expanding the so- 
called quantum determinant of the transfer matrix [45] . 

As mentioned before, the energy level statistics of Haldane-Shastry type 
models is neither Poissonian nor of Wigner type. We propose to understand 
this fact via the observation that this family of models fall in the class of quasi- 
polynomial quantum integrable, rather than linear or polynomial. 

Richardson-Gaudin type models 

Yet another family of models which are interesting to discuss are the Richardson- 
like models [53] [SH [55] [56] [57]. Similarly to the Haldane-Shastry- type mod- 
els, the Hamiltonian has quadratic density character. However, the conserved 
charges are simpler [SB], and are of linear density character. We can thus put 
this family of models in the class of linear quantum integrable in the sites basis. 

Finally, in a generic model (which still fulfills QI:N), it can be expected 
that all charges will be of exponential density character. We would thus by con- 
vention say that if a model is not at least sub-exponential integrable, it simply 
should be called non-integrable. 

We could carry on mentioning other specific models. However, for the sake 
of compactness, and given the fact that the examples we have given clearly 
illustrate our idea, we now turn to the more crucial question of the physical 
interpretation of the new definition. 

6 Physical meaning of the new definition 

One of our requirements was that the classification scheme should be physically 
meaningful. We have already mentioned a few features traditionally associated 
with integrable models, such as factorized scattering, nondifFraction, particular 
energy level statistics, etc. The main theme which we would however like to 
understand better is how the classical notion of (non)ergodicity translates to 
the quantum case. 

In this context, the appropriate starting point for discussion is probably 
Mazur's inequality [58], which has already proven useful in a number of cases, 
for example in [SI] . Let us recall its statement [SHI (SSI ISI] • Consider an operator 
A, and consider the canonical average (Ayl(O)AA(t)) where AA = A — {A). If 
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{Sq}, oi — 1,...,Nq is a set of constants of motion, the long-time average 
correlation obeys 

.'if. f *{AA(0) A^«» > I !«^^ (9) 

in which the charges are defined in such a way as to be 'orthogonal' for the 
chosen form of averaging, 

{QcQp) ^ S,,p{Ql) (10) 

(note that we write the charges here as Qa to distinguish them from the charges 
we have used in Section|31 our charges do not depend on the averaging procedure 
chosen, while the Qa charges do by implementation of (ITUl) '). The inequality 
becomes an equality if the set of charges used is maximal. The operator A is 
called ergodic if 

lim 1 / dt{AAiO)AA{t)) = 0, (11) 

in other words if it is orthogonal to all the conserved charges. The physical 
picture is that all correlations have dissipated in Hilbert space, and no quasi- 
periodic behaviour can be identified. Note the fact that in order to prove non- 
ergodicity, it is sufficient to find a single conserved charge having a finite overlap 
with the operator of interest (more elaborately, an infinite set of vanishing con- 
tributions summing up to a finite value) , thereby making the right-hand side of 
([9]) nonvanishing. 

How does our definition of integrability relate to Mazur's inequality? Let us 
make the following reasoning. We can first of all generalize Mazur's inequality 
to a generic energy-diagonal expectation value, 

(...) ^ (...)f ^ ^5]/,(a|...|a), ^ (12) 

a a 

in which a labels eigenstates and fa are specified real parameters. The form 
of the distribution is important: Mazur's inequality trivializes to an empty 
statement if the average is taken over a single state, since the conserved charge 
eigenvalues then cancel in the numerator and denominator in each term of the 
right-hand side. We thus focus our attention on distributions with at least some 
measure of extensivity. 

Let us consider physically meaningful operators A which are represented by 
simple matrices in the preferred basis. Such operators would thus be represented 
as 

A,,...,„^e^--"" (13) 



with Ne {A) (using the terminology of Section |4]) being for example such that A 
is of polynomial density character (this would thus include all operators made 
of finite products of on-site operators or the Fourier transform of such objects). 
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Let us now consider computing the building blocks of Mazur's inequality, 
namely the expectation values 

((AA)Q„)f. (14) 

This is really a measure of the 'overlap' of the operator AA with the conserved 
charges we choose to work with. Following (SD], this can be translated at the 
operator level for an 'orthogonal' (obeying ([TU])) set of charges Qa- One can 
write 

A = (15) 

a 

where Oq, are c-numbers (implicitly depending on the averaging scheme chosen) 
and A' is an operator which is off-diagonal in energy, and therefore doesn't 
contribute to Mazur's inequality. The ratios in ^ are then simply related to 
the coefficients 

((AA)g„)f 

The next steps in the discussion depend on the integrability class our system 
falls into. Let us discuss the generic non-integrable case first. In this case, 
all the conserved charges Qa are expressed in the preferred basis in terms of 
dense matrices with 0(dim(?^)^) entries, and have thus exponential density 
character. Pick any one of these charges, say Qai, and use it as the basis for 
constructing the set Qa of orthogonal charges to be used. The exponential 
density character of Qa^ charge means that the coefficient Ua^ in (fTSj) will be 
exponentially small. Picking a second charge Qq2: using an orthogonalization 
procedure to satisfy (fTOj). one remains (except for utterly contrived cases) with 
an exponential density character charge Qa2 having vanishing overlap with 
our operator A. This simply carries on, and no charge can be found which 
overlaps with A, so Mazur's inequality indicates ergodicity. 

Let us now turn to the case of an integrable model according to our defini- 
tion (for example, a linear quantum integrable model, for which we have many 
charges Qa at our disposal which can be written as single summations of local 
operator products). By their nature, these charges can be expected to have a 
finite overlap with operator A. Enforcing the orthogonality condition ([T0| does 
not change the density character of the charges, so we obtain a set of orthogonal 
charges Qa whose (summed) overlaps have the potential to give a finite value 
to the right-hand side of Mazur's inequality. This must of course be checked 
in individual cases, but the door is clearly left open by the conserved charges' 
simple structure, related (in the sense of their simple density character) to that 
of operator A. The same conclusions would hold for other integrability classes, 
for appropriate operators having a sufficient degree of similarity to the (simple) 
conserved charges. 

Let us make some simple and obvious remarks. First of all, of course, er- 
godicity only makes sense in the infinite system size limit. For a finite system, 
in parallel to the naive definition QI:N of quantum integrability, any quantum 
system will possess 'reasonable' observables failing to show ergodicity in Mazur's 
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sense; these observables will however not be physically meaningful for what we 
have referred to as non-integrable systems, since they must have exponential 
density character. One thus cannot say that a system is ergodic or not, but 
rather only that a given observable A is, for a given averaging in a given sys- 
tem. Our definition offers a way of understanding the origin of this ergodicity, or 
of its absence, by looking at the detailed structure of the conserved charges, and 
their relationship with the observables considered within the adopted averaging 
scheme. Our scheme of separating models into different integrability classes in 
fact leads us to expect that ergodicity manifests itself differently in each of these 
classes. We will further characterize this in subsequent publications. 

Another context in which the role of quantum integrability is frequently 
discussed is the dynamics resulting from a non-equiblibrium initial state [62l 
I63[ I64i I65j . A common way to bring a system out of equilibrium is by means of 
a so-called quantum quench j66j . Consider a Hamiltonian H\ which depends on 
some parameter A. The quantum quench now consists of preparing the system 
in an eigenstate of Hx, followed by instantaneously changing the parameter 
A — ?• A'. If [H\,Hx>] ^ 0, one expects complicated non-equilibrium dynamics 
to ensue. To clarify this, let \^) be an eigenstate of H\ before the quench. 
We can now express |(/>(i)) for i > as e~'*^^'|(/)) = X^n "^"1")^"'^"*' where 
the summation is over all eigenstates of and the coefficients c„ — {n\(\)) 
represent the overlaps between the initial state |(^) and the eigenstates \n) of 

. The complicated form of \(\>{t)) makes it hard to make predictions at short 
time scales. It is common to simplify the problem by considering the long time 
average of an observable A, 

hm i C m\A\m) = E VnWn\A\n) (17) 

^' n 

which has the same energy-diagonal form as (|9])Q For most quenches, we can 
also make the assumption that the distribution {|c„p} is narrow in energy. For 
a sufficiently large system one can try to approximate the distribution {|cnp} 
by a canonical ensemble 

Y,\c^\^n\A\nf ^^TriAe-^"} (18) 

n 

where the inverse temperature P is determined from the initial condition 

{m\H\m) = ^'^^{He-f'"] (19) 

and Z is the usual partition function Z = Trje"*^^}. This approximation, while 
successful for some cases, fails in general for integrable systems, this failure being 
typically attributed to the presence of additional conserved charges. It is thus 

*As in Mazur's inequality, in general the summation also contains ofT-diagonal expectation 
values (n|A|m,) as a result of degenerate energy levels. We neglect these here for simplicity. 
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natural and appealing to generalize the canonical ensemble by the so-called 
Generalized Gibbs ensemble [6 7) , 

{A)GGE = ^Tr{Ae-^^^^^^}, (20) 

the generalized partition function Z being now computed as Z — Tre~^>^'^'. 
Initial (quench-time) conditions fix the Lagrange multipliers via the self- 
consistency requirements 

immm) = ^Tr{Q,e-^.^^2^}. (21) 

Using all available conserved charges as per the original formulation of this 
approach (in other words by adopting the QI:N view of integrability) , one might 
consider the eigenstate projectors Q„ = |?T-)(fi| as a maximal set of conserved 
charges to be exploited. While this would lead the an exact description, this is of 
limited use since in this case e~^"^" = |c„|^ are simply the overlap coefficients 
mentioned above, this involving 0{dim{H)) parameters. What lacks in the 
description of the Generalized Gibbs ensemble is thus a prescription for which 
conserved charges should actually be included in the scheme. Since we assumed 
that the distribution {|c„p} is narrow in energy, one should include (besides 
the inevitable Hamiltonian) conserved charges that have a small variance in the 
distribution {|c„p}. In contrast, if a conserved charge Q„ has large variance, 
its Lagrange multiplier /?„ would be negligible, and including this charge in the 
scheme would not provide meaningful state selection. For an integrable system 
as defined in Section 4, many conserved charges having a similar character as the 
Hamiltonian are available (these having non-negligible Lagrange multipliers), 
and one cannot generally neglect their effect on the distribution. On the other 
hand in the non-integrable case most of the conserved charges will have an 
exponentially large variance for this particular distribution, and can thus be 
disregarded. Our proposed definition of quantum integrability thus leads us to 
suggest using restricted generalized Gibbs ensembles, which we will return to in 
future work. 

7 Conclusions 

One of the difficulties in discussing quantum integrability is of course that it 
can very quickly feel at best confusing, and at worst rather pedantic. We feel 
that these Proceedings offered the right kind of forum for this contribution, in 
view of the importance of an increased clarity of concepts for use in discussions 
about quantum many-body equilibration dynamics. Instead of attempting to 
define 'quantum integrability', it might have been more reasonable and less 
confusion-inducing to simply introduce new terminology avoiding the keyword 
'integrability' altogether. We have tried to mitigate this risk by suggesting 
to associate the terminology of algorithmic time complexity with any specific 
mention of quantum integrability. If we say then that the Heisenberg chain is 
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linear quantum intcgrablc, it will be clear that the notions and definition of our 
paper are the ones that are referred to, and which consequences we associate to 
such a categorization. 

One point which wc should clarify is how wc reconcile the current concepts 
with Integrable Field Theories (ITF) and Conformal Field Theories. The point 
is that we have considered in detail how the morphology of the conserved charges 
of a given model behaves as the system size is increased, i.e. as one moves to- 
wards a continuum/field theory limit. If a model is classified as linear integrable 
in the sites basis within our scheme, it will automatically lead to a field theory 
having an infinite number of conserved charges expressible as integrals of local 
densities by choosing an o{N) subset of charges in the finite lattice regulariza- 
tion. 

In the beginning, we listed our three main requirements: that a proper 
definition of quantum integrability should be unambiguous, and should lead to 
different classes of models, these having different discernible types of physical 
behaviour. On this last point, while we have provided heuristic arguments, the 
specific and detailed consequences of our definition of quantum integrability 
have not yet been fully worked out to the level of general theorems. This first 
involves a rather immense amount of work of looking at specific examples in 
detail, which we are pursuing at this time. This will teach us whether our 
definition actually has interesting substance, which is perhaps not completely 
convincing at this point. While we think that our definition is not ambiguous, 
the level of mathematical rigour we have adopted is rather low in comparison 
to what is achievable in the field of integrable models, but should be put in 
correspondence to the level of discussion on algorithmic time complexity. We can 
thus be satisfied overall that our initial objectives feel as if they have been met. 
Much work remains to be done, and we will extract more specific conclusions 
on particular classes of models and observables in future publications. 
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